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Abstract 

The first moments of the polarized valence parton distribution functions (PDFs) trun- 
cated to the wide Bjorken x region 0.004 < x < 0.7 are directly (without any fitting 
procedure) extracted in the next to leading order (NLO) QCD from both COMPASS and 
HERMES data on pion production in polarized semi-inclusive DIS (SIDIS) experiments. 
The COMPASS and HERMES data are combined in two ways and two scenarios for the 
fragmentation functions (FFs) are considered. Two procedures are proposed for an es- 
timation of light sea quark contributions to the proton spin. Both of them lead to the 
conclusion that these contributions are compatible with zero within the errors. 
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Longitudinally polarized parton distribution functions (PDFs), and especially their first 



o ■ moments, which directly compose the nucleon spin together with the orbital parton momenta, 
p , are of crucial importance for solution of the proton spin puzzle, attracting great both theoretical 
<D \ and experimental efforts during many years. Nowadays, there is a huge growth of interest to 
£j| the semi-inclusive DIS (SIDIS) experiments with longitudinally polarized beam and target such 
as SMC m, HERMES 0, COMPASS (see, for instance, g] for review). It is of importance 
j> ■ that the SIDIS experiments, where one identifies the hadron in the final state, provide us with 
^ . the additional information on the partonic spin structure in comparison with the usual DIS 
^ ! experiments. Namely, in contrast to the DIS data, the SIDIS data allows us to find the sea 
"xl" \ and valence contributions to the nucleon spin in separation. In this paper we just focus on this 
' important task. To this end we apply the original procedure of the polarized SIDIS data analysis 
^ ■ in NLO QCD, elaborated in the sequel of papers [SJ El CTJ (see also [I] for more details). It is 
t-h . of importance (especially for analysis of still relatively poor SIDIS data we deal with) that this 
J> ! alternative procedure allows to extract the truncated to the accessible for measurement Bjorken 
\ x region moments of polarized PDFs directly, without any fitting procedure (with unavoidable 
5^ ■ arbitrariness in the choice of functional form of Aq at initial Q 2 and a lot of free varying 
parameters): within this procedure the central values of asymmetries and their uncertainties 
directly propagate to the extracted values of moments and their errors. Of importance also 
that we use in our analysis the difference asymmetries ^just as in Refs. [5j [7]), that allows to 
avoid the application of poorly known FFs (such as Dff and Dg FFs) - see [SJ [7] and also jl] 
for review. 

Within this paper we use two NLO parametrizations AKK08 |8J and DSS [H] for FFs, which 
differ quite essentially in the pion sector. AKK08 parametrization corresponds to unbroken 
577/(2) symmetry: 

D x = Dl + £ Df SU ^ 2) Df £ Df, (1) 
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D 2 = Df = Df SU ± 2) Df = Df . (2) 

On the other hand, DSS parametrization [U] allows violation of SUf(2) symmetry in the sector 
of the favored pion FFs: 

D 1 = Df^Df, (3) 
Di = Df = Df, (4) 

D2 = Df = Df 5 = 2 ^ Df = Df , (5) 

so that the favored FFs D\ and D\ of u and d quarks are not equal to each other. 

For the scenario preserving SUf(2) symmetry (AKK08 parametrization for FFs here) the 
procedure of direct extraction in NLO QCD of n-th moments of the valence PDFs from the 
measured difference asymmetries is described in Refs. [5], [7] in details. Let us recall the key 
necessary equations. 

The theoretical expressions^] for the difference asymmetries in NLO QCD look as 



Ap n (x, Q 



2-. 



z f z dz h {Au v -d v )k{D 1 -D 2 ) ' 



for the proton target, and 



= (i + R) (i - L n ) J * d T {Aur + Ad ^ Ak ^ ~ D ^ , (7 ) 

z 2 Sldz h (u v + d v )K(D l -D 2 ) 

for the deuteron target. Here, for brevity, we have introduced the operator notation K and 
AK: 



K = l + ®^C 2 qq ®, AK = 1 + ®^AC qq ®, 

2TT ATI 

so that 



[ q KD)(x,z h )= q (x)D(z h ) + ^ r^^^^d), 0) 

J X J Zfr 



[A q AKD](x,z h ) = A q (x)D(z h ) + ^- f % f %Aq (-) AC qq {x' , z')D (^) . (10) 

Z7T J x X J Zh Z \X / \ Z / 

In Eqs. ()6]), (JU) (1 — 1.5ud), ojd = 0.05 ± 0.01, is the factor accounting for the deuteron 
D-state contribution, the quantity R = <7l/<Jt is taken from Ref. [12], AC qq is the Wilson 
coefficient entering NLO QCD expression for the polarized SIDIS structure function g\ and can 
be found, for instance, in [T3], C q = C qq + C qq is the Wilson coefficient entering the NLO QCD 
expression for the unpolarized SIDIS structure function F^, while the coefficients C qq and C qq 
(entering F± and F^) also can be found in [13J. In our subsequent calculations we also use 
NLO parametrization GJR08 [H] for unpolarized PDFs. 

4 Since namely F% function is measured in experiment and then is used to parameterize unpolarized PDFs, 
the form of these equations used here is the most convenient to properly account for correction due to the factor 
R = ctl/pt (see discussion around Eq. (10) in [10] , around Eq. (12) in [TTJ and references therein). 
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The equations allowing to find from the data on difference asymmetries the n-th moments 
A n q = f a dx x n ~ 1 q(x) of valence PDFs truncated to the accessible for measurement x region 
[a, b] look as 

A n U V ~--P -i-; A n dy~--^ (11) 

MnH — Mn)2 Mti)l — Mn)2 



where 



.4 



(n) 



(1- I.Swd)- 1 



ttear n - 1 (l + R)' 1 / d2 A [(4« v - d v )K(D l - D 2 )], (12) 



dxx n -\l + R)- 1 dz h [{u v + d v )K{D l -D 2 )}. (13) 



The quantities L( n )i, £( n )2 are defined as 



■ 7T — 

J (n)d 



£(n)l 


r 7T T 

= ^(n 


£(n)2 


r tt^ 

= L {n 




-/ 

Jz 



' 7V^~ _ 7" 7T 



J (n)« 



J (n)u> 



J z h 



where 



/ dx x n 1 AC qq (x, z). 
Jo 



(14) 



A lot of numerical tests performed in Refs. [7] have shown that Eqs. (fTTj) allow to 
reconstruct the truncated moments of polarized valence PDFs with a high precision even in the 
case of rather narrow HERMES x region 0.023 < x < 0.6, divided on relatively small number 
of bins (nine bins only) . 

It is easy to see that for scenario with the broken SUf(2) symmetry one has instead of dSJ), 
(J7|) the equations 



A TT — 7T 



{AAu v - Ad v )AK{D 1 - D 2 ) + Ad v AK{D l - D 1 
{Au v - d v )K(D 1 - D 2 ) + d v K(D 1 - Dx) 



^ PV1 3 (A^y + Ad v )AK((D 1 - D 2 ) + |px - gQ) 

1 + -fill 1 Wn ) = — r 

2 ( Wy + dy)K((A-£>2) + f (£>!-£>!)) 



(15) 

(16) 



The generalization of Eqs. (fIT]) - (fl4"]) on the case of broken SUf{2) symmetry (Eqs. (J!])-®) is 
also straightforward. 

Both COMPASS USDS] and HERMES [2] collaborations published the data only on asym- 
metries Al d , while the published data on the pion difference asymmetries Ap d ~ n are stilfl 



3 At present this work is in preparation at COMPASS. 
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absent. That is why the special procedure was applied in [7] to construct asymmetries A^ -77 
from the HERMES data on pion production, and we repeat here this procedure for the COM- 
PASS case. Namely, in each z-th bin the pion difference asymmetries can be rewritten as 

A« + -*-{ Xi ) = fj' A* + ( Xl ) - A*-( Xi ), (17) 

R { - 1 R t - 1 

where R% is the ratio of unpolarized cross-sections for 7r + and tt~ production: R^~ = 
a unpoi( x O/ a un P oi( x i) = / N? . As it was argued in Ref. [7] this relative quantity is very 
well reproduced by the the LEPTO generator of unpolarized events [TTJ, which gives a good 
description of the fragmentation processes. So, we again use here the LEPTO generator for 
this purpose. 

Let us now discuss the question of Q 2 dependence of asymmetries and its influence on the 
final results. The point is that both DIS and SIDIS asymmetries very weakly depend on Q 2 
(see, for instance, Fig. 5 in Ref. [IS]), so that the approximation 

A(x u Q 2 )~A(x u Q 2 ) (18) 

is commonly used (see, for example, Refs. [H El EES]) f° r analysis of the DIS and SIDIS asymme- 
tries. Nevertheless, for more comprehensive analysis, it is useful to account for the corrections 
caused by the weak Q 2 dependence of the difference asymmetries, i.e., to estimate the shifts 

5 t A;^- = a£t~ ( Xii Ql) - a;^- ( Xi , q 2 ) (19) 

in the difference asymmetries and their influence on the moments of valence PDFs. To this 
end we first approximate r.h.s of Eq. (|T9|) by the respective difference of "theoretical" asym- 
metries calculated with substitution of two novel parametrizations [THJ, [TT] on polarized PDFs 
(elaborated with application of both DIS and SIDIS data) into the NLO QCD equations fl§D, 
flTJ (or f|ToT) - ffToT) ). and then averag^l the obtained results on SiA^f 77 . Adding the calcu- 
lated in this way 5iAp + d ~ w to the initial experimental asymmetries A 77 ^ 77 (xi,Q 2 ), we esti- 
mate the evolved from Q 2 to Q\ asymmetries A 77 ^ 77 (x{, Q%)\ e voived- Using the obtained in 
such a wajQ evolved asymmetries we extract the respective corrected moments of the valence 
PDFs A' n q v I corrected- Then we compare the corrected moments A' n q v \ corrected with the respec- 
tive moments A' n q v obtained within the approximation (IT51) . and calculate the respective shifts 
S(A' n q v ) = A' n q v \ corrected — A! n q v as well as the relative quantities 5(A' n q v ) / A' n q v . 

Of importance is the optimal choice of the scale Ql common for evolved asymmetries, 
allowing as much as possible to reduce the shifts in results due to evolution. Our experience 
shows that for the combined analysis of COMPASS and HERMES data (see below) the optimal 
choice is close to Ql = 10 GeV 2 . 

The evolved to Q 2 = 10 GeV 2 pion difference asymmetries constructed with Eq. (1171) from 
the COMPASS [H Q5] and HERMES [2] data on A^ d SIDIS asymmetries are presented in 
Figs. \1\ and [2] respectively. 



6 Notice that the shifts in asymmetries as well as in the final results on the moments of valence PDFs obtained 
with two applied parametrizations differ very insignificantly from each other. 

7 Notice that the considered procedure of the asymmetry evolution is quite similar to the procedure used by 
SMC for the ^i P (d) reconstruction (see Section YB in Ref. [20]). Notice also that the original procedure of Ax 
evolution was elaborated in Refs. [2T1I22] . 
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We perform the combined analysis of COMPASS [TBI [15] and HERMES [2] data on pion 
production with both proton and deuteron targets. COMPASS collaboration published their 
data on A* d in the Bjorken x ranges 0.004 < x < 0.7 and 0.004 < x < 0.3 for the proton and 
deuteron targets, respectively, while the HERMES data on these asymmetries were presented in 
the range 0.023 < x < 0.6 for both targets. The statistical summation of asymmetries Ap d ~ n 
(constructed with (II 7ft) is performed in accordance with the standard formulas 



Ah I A%\ exp i/ (5A%\ exp i) 2 + A%\ exp 2/(5Apf \exp2) 2 / 9f .N 



(SA N \ averaged ) - 1/{5A h nUpi)2 + i/ ( ^| e ^ 2) 2- ( 21 ) 

At the same time one can apply Eqs. (1201 . (12TI) directly only for coinciding x bins of different 
experiments. However, this is the case only for last three bins of COMPASS and HERMES 
experiments we deal with (after proper extrapolation^] of HERMES data in the last bin from 
0.6 to 0.7 upper x value). Besides, notice that for two last bins the COMPASS published SIDIS 
data for deuteron target are still absent. That is why it is of especial importance (and we do 
it first of all) to include in the analysis of COMPASS data the HERMES data in the region 
0.2 < x < 0.6 (last three bins of HERMES). The results on the difference asymmetries obtained 
in such a way are presented in Fig. [3J 

The respective results on the moments of polarized valence PDFs are presented in Table [TJ 



Table 1: Four first moments of polarized valence PDFs truncated to the region 0.004 < x < 0.7 are presented 
at Q 2 = 10 GeV 2 . The moments are obtained as a result of NLO QCD analysis of the combined data on A* d _7r 
(see Fig. constructed with Eq. JT^I from the COMPASS data on A^ d in the regions 0.004 < x < 0.7 (proton 

target), 0.004 < x < 0.3, (deuteron target), and HERMES data on A* d in the region 0.2 < x < 0.6 (last three 
bins of HERMES). Capital letters A and B correspond to the application of AKK08 and DSS parametrizations 
for FFs, respectively. Rome numbers I and II correspond to the moments uncorrected and corrected due to 
evolution, respectively. Besides, the relative corrections 5 r (A' n qv) = 8(A' n qy)/ A' n qy caused by evolution are 
presented here. 



A' n u v 


n 


A: 


A n 


5 r (A' n u v ) 


Bi 


Bn 


5 r (A' n u v ) 


1 


0.731 ± 0.087 


0.695 ±0.087 


-3.8% 


0.693 ±0.084 


0.713 ± 0.084 


2.8% 


2 


0.166 ±0.024 


0.167 ±0.024 


0.8% 


0.155 ±0.024 


0.158 ±0.024 


1.6% 


3 


0.055 ±0.010 


0.055 ±0.010 


1.3% 


0.052 ±0.010 


0.052 ±0.010 


1.8% 


4 


0.022 ±0.005 


0.022 ±0.005 


1.5% 


0.021 ±0.005 


0.021 ±0.005 


2.0% 


A' n d v 


n 


A: 


A n 


Sr(Kd V ) 


Bi 


B n 


5 r {A'J v ) 


1 


-0.519 ±0.162 


-0.524 ±0.162 


0.9% 


-0.473 ±0.157 


-0.481 ±0.157 


1.7% 


2 


-0.100 ±0.054 


-0.102 ±0.054 


1.8% 


-0.090 ±0.051 


-0.092 ±0.051 


2.7% 


3 


-0.029 ±0.023 


-0.030 ±0.023 


2.5% 


-0.026 ± 0.022 


-0.027 ±0.022 


3.7% 


4 


-0.011 ±0.011 


-0.011 ±0.011 


3.1% 


-0.010 ±0.010 


-0.010 ±0.010 


4.4% 



8 Our experience show that this leads to negligible changes in the final results, irrespective of the choice of 
the extrapolation procedure. So, we apply here the simplest way of extrapolation, prescribing to A^ d | H ERMES 
to be in the extended region [0.4, 0.7] the same as in the last HERMES x-bin [0.4, 0.6]. 
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On the other hand, to maximally increase the available statistics it is, certainly, very de- 
sirable to perform the complete combined analysis, using the COMPASS and HERMES data 
taken in the entire x regions accessible for these experiments. To this end we have elaborated 
the special procedure allowing to combine the data on asymmetries coming from experiments 
with different binnings - see the Appendix. The results on the difference asymmetries obtained 
in such a way are presented in Fig. HJ 

The respective results on the moments of polarized valence PDFs are listed in Table [2J 

Table 2: Four first moments of polarized valence PDFs truncated to the region 0.004 < x < 0.7 are presented 
at Q 2 = lOGeV 2 . The moments are obtained as a result of NLO QCD analysis of the combined data on 
A p*d"*~ ( see Fi §- BI, constructed with Eq. (JT7J) from the COMPASS and HERMES data on A^ d in the entire 
x-regions accessible for measurement (0.004 < x < 0.7, 0.004 < x < 0.3 for COMPASS and 0.023 < x < 0.6 for 
HERMES). Capital letters A and B correspond to the application of AKK08 and DSS parametrizations for FFs, 
respectively. Rome numbers I and II correspond to the moments uncorrected and corrected due to evolution, 
respectively. Besides, the relative corrections 5 r (A' n qv) = 5(A' n qv)/ A' n qv caused by evolution are presented 
here. 



A' n u v 


n 


A: 


A n 


5 r (A' n u v ) 


Bi 


B„ 


S r (A' n u v ) 


1 


0.712 ± 0.078 


0.660 ±0.078 


-4.3% 


0.683 ±0.076 


0.711 ±0.076 


4.0% 


2 


0.166 ±0.023 


0.168 ±0.023 


0.8% 


0.156 ±0.024 


0.159 ±0.024 


1.7% 


3 


0.055 ±0.010 


0.056 ±0.010 


1.4% 


0.052 ±0.010 


0.053 ±0.010 


1.8% 


4 


0.022 ± 0.005 


0.022 ±0.005 


1.6% 


0.021 ±0.005 


0.021 ±0.005 


2.0% 


A' n d v 


n 


Ai 


A n 


Sr(Kdv) 


Bi 


B n 


5 r {A'J v ) 


1 


-0.414 ±0.149 


-0.427 ±0.149 


3.0% 


-0.376 ±0.145 


-0.381 ±0.145 


1.4% 


2 


-0.087 ±0.053 


-0.089 ±0.053 


2.5% 


-0.078 ±0.051 


-0.080 ±0.051 


2.7% 


3 


-0.027 ±0.023 


-0.028 ±0.023 


2.8% 


-0.024 ±0.022 


-0.025 ±0.022 


3.8% 


4 


-0.010 ±0.011 


-0.011 ±0.011 


3.3% 


-0.009 ±0.010 


-0.010 ±0.010 


4.6% 



In Tables [HE] both scenarios Eqs. (HJ, (AKK08 parametrization) and Eqs. ©-© (DSS 
parametrization) for FFs are considered. One can see that the results are consistent within the 
errors. Besides, for each scenario we present the results obtained with and without corrections 
due to evolution of asymmetries. It is seen that the difference is not too significant (the relative 
corrections 5(A' n q v ) / A' n q v take the small values). 

Thus, we estimated in NLO QCD the contributions of valence quarks (first moments of 
polarized valence PDFs) to the nucleon spin. Let us now estimate the respective contributions 
of light sea quarks. We do it in two different ways. 

Within the first procedure one first of all uses some NLO QCD parametrization on the 
polarized PDFs to estimate the quantities A[q + A[q (q = u,d). Since the sums Aq(x) + 
Aq(x) (q = u, d) are well fitted by the precise purely inclusive DIS data (these quantities 
are considered as relatively well known and practically are the same for the different modern 
parametrizations) it is not especially important which parametrization one applies for this 
purpose (here we use the most popular and widely cited DSSV (19] parametrization). Then, 
having both (A[q + A^g)| parametrization (q = u,d) and (see Tables Q] and |2]) A[q v (q = u,d) 
quantities one easily gets the truncated first moments of sea u and d quarks, applying the 
obvious relation 

A' l q= l -[(A' l q + A' l q)-A' 1 q v }. (22) 



6 



The first moments A[u, A[d obtained in this way, as well as their differences and sums are 
presented in Tables [3] and HJ 



Table 3: First moments of polarized sea PDFs truncated to the region 0.004 < x < 0.7 are presented at 
Q 2 = lOGeV 2 , as well as their sums and differences. The moments are obtained with application of Eq. (|2"2")l . 
where DSSV parametrization is used to estimate (A^g + A'^lparamctrization, while the first moments of valence 
PDFs are taken from Table [T] (HERMES data only from three last bins are applied) . Capital letters A and B 
correspond to the application of AKK08 and DSS parametrizations for FFs, respectively. Rome numbers I and 
II correspond to the moments uncorrected and corrected due to evolution, respectively. 





Ai 


An 


Bi 


B n 


A[u 

A[d 
A[u + Aid 
A[u - Aid 


0.018 ± 0.044 
0.065 ±0.081 
0.082 ±0.092 
-0.047 ±0.092 


0.036 ± 0.044 
0.067 ±0.081 
0.102 ±0.092 
-0.032 ±0.092 


0.037 ±0.042 
0.042 ±0.079 
0.078 ±0.089 
-0.005 ±0.089 


0.027 ±0.042 
0.046 ± 0.079 
0.072 ±0.089 
-0.019 ±0.089 



Table 4: First moments of polarized sea PDFs truncated to the region 0.004 < x < 0.7 are presented at 
Q 2 = 10 GeV 2 7 as well as their sums and differences. The moments are obtained with application of Eq. ((22)) . 
where DSSV parametrization is used to estimate (A' 1 q + A'^lparametrization, while the first moments of valence 
PDFs are taken from Table [5] (all data on pion production of COMPASS and HERMES are combined). Capital 
letters A and B correspond to the application of AKK08 and DSS parametrizations for FFs, respectively. Rome 
numbers I and II correspond to the moments uncorrected and corrected due to evolution, respectively. 





Aj 


A n 


Bi 


Bn 


A[u 

A[d 
A[u + A[d 
A[u - A[d 


0.027 ±0.039 
0.012 ±0.075 
0.039 ±0.084 
0.015 ±0.084 


0.053 ±0.039 
0.019 ±0.075 
0.072 ± 0.084 
0.034 ±0.084 


0.042 ± 0.038 
-0.007 ±0.073 
0.035 ±0.082 
0.048 ± 0.082 


0.028 ±0.038 
-0.004 ±0.073 
0.023 ±0.082 
0.032 ±0.082 



The idea of alternative procedure for investigation of the sea contributions to the nucleon 
spin in NLO QCD is based on the proper application of SUf{2) (Bjorken sum rule) and SUf(3) 
sum rules: 



«3 



a 8 



(Aiu + Aim) - (A^ + A^) 



9a 



9v 



F + D = 1.2670 ±0.0035, 



(23) 



Aiw + Aiw ± Aid + Aid - 2(A 1 s + A ± s) = 3F — D = 0.585 ± 0.025. (24) 



Bjorken sum rule (l2"3"j) rewritten in terms of valence and sea distributions produces quite 
$ood approximation [5j [23] (see Ref. jl] for review) 



Aiw 



9a 
9v 



~ (A[u v - A[d v ) 



(25) 



for the difference of full (not truncated!) moments Aiu and Aid even in the case of rather narrow 
HERMES x-range, while for the wide COMPASS x-range we deal with here this approximation 
works very well - see the respective numerical tests in Ref. [3]. The point is that since the 
valence PDFs (contrary to the sea PDFs) are suppressed near the low boundary x = (so that 
Ag = Aqv + Aq — > Aq as x — > 0), the omitted in r.h.s. of Eq. (125]) term f^ dx (Auy — Ady) is 
small even for HERMES low x boundary a = 0.023, and becomes really negligible for COMPASS 
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a = 0.004 we deal with here. In turn, another omitted term J Q7 dx(Au\r — Ady) in r.h.s of 
(I2"5j) is also negligible since all PDFs just die out at so high x values. 

On the other hand, to estimate the sum of full moments in NLO QCD we use the sum 
rule ()24p and purely inclusive DIS data on the first moment Tf of deuteron structure function 
gid (measured with high precision). To this end we use the NLO QCD expression for = 

rf ee (1 - l.Su^-M = \{V\ + r?) = (l - 2^1) Q- 6 a 8 + ^(Q 2 )) , (26) 
which produces very good approximation for AiU + Aid: 

Aiu + Aid ~ ^3(l + ^)rf ±^s) -i(A> y + Aid y ), (27) 

where we again omitted the small contributions of valence PDFs / °' 004 dx (Auy + Ady) and 
f Q7 dx (Auy + Ady)). We use in Eq. f )27|) the numerical value of rf taken from the COMPASS 
paper [TO] : 

Tf = 0.051 ±0.003 ±0.006. 

The obtained with Eqs. (125]) . (12 7p results on the sums and differences of the first moments 
of sea PDFs, as well as on the moments themselves in separation are presented in Tables Eland 
El 

Table 5: Sums and differences of the first moments of polarized sea PDFs, as well as the moments themselves, 
obtained in NLO QCD at Q 2 — 10 GeV 2 within the approximations f|25[) and (127)) . The truncated first moments 
of valence PDFs are taken from the Table [T] (HERMES data only from three last bins are applied). Capital 
letters A and B correspond to the application of AKK08 and DSS parametrizations for FFs, respectively. Rome 
numbers I and II correspond to the moments uncorrected and corrected due to evolution, respectively. 





A: 


An 


Bi 


Bn 


Aiu 

Aid 
Aim + Aid 
Aiu - Aid 


0.059 ± 0.045 
0.050 ±0.082 
0.109 ±0.095 
0.009 ±0.092 


0.077 ± 0.045 
0.053 ±0.082 
0.129 ±0.095 
0.024 ± 0.092 


0.078 ± 0.043 
0.027 ±0.079 
0.105 ±0.092 
0.051 ±0.089 


0.068 ± 0.043 
0.031 ±0.079 
0.099 ±0.092 
0.037 ±0.089 



Table 6: Sums and differences of the first moments of polarized sea PDFs, as well as the moments themselves, 
obtained in NLO QCD at Q 2 = 10 GeV 2 within the approximations (|25p and (f!?T)) . The truncated first moments 
of valence PDFs are taken from the Table [5] (all data on pion production of COMPASS and HERMES are 
combined). Capital letters A and B correspond to the application of AKK08 and DSS parametrizations for FFs, 
respectively. Rome numbers I and II correspond to the moments uncorrected and corrected due to evolution, 
respectively. 





Ai 


An 


Bi 


B n 


Aiu 

Aid 
Aiu + Aid 
Aiu - Aid 


0.068 ±0.041 
-0.002 ±0.075 
0.066 ±0.087 
0.071 ±0.084 


0.094 ±0.041 
0.004 ±0.075 
0.099 ±0.087 
0.090 ±0.084 


0.083 ± 0.040 
-0.021 ±0.073 
0.061 ±0.085 
0.104 ±0.082 


0.069 ± 0.040 
-0.019 ±0.073 
0.050 ±0.085 
0.087 ±0.082 
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Looking at Tables El H and Tabled U E one can draw an unexpected conclusion, that 
irrespective of the procedure used in the SIDIS data analysis the first moments of sea PDFs 
are consistent with zero within the errors. In particular, in contrast with the different model 
predictions [2JJ [25] |26j [27] EH1 EH] the polarized sea asymmetry Aiu — Aid appear to be just 
zero within the errors. It is of importance, because some of these models (see Refs. [25] [27] 
and references therein) predict that this asymmetry is even larger in absolute value than the 
unpolarized sea asymmetry u — d (which takes quite considerable value allowing to explain a 
large violation of the Gottfried sum rule). 

In conclusion, let us briefly discuss the obtained results. 

The first moments of the polarized valence PDFs truncated to the wide Bjorken x region 
0.004 < x < 0.7 are directly (without any fitting procedure) extracted in NLO QCD from both 
COMPASS and HERMES polarized SIDIS data. To this end we apply two scenarios for the 
fragmentation functions and two ways to combine the COMPASS and HERMES data on pion 
production (one of them is based on the proposed special procedure, allowing to combine all 
data coming from experiments with different binnings). In turn, the obtained results on the 
valence PDFs have allowed us to estimate in two ways the contributions of light sea quarks to the 
proton spin which, surprisingly, occur compatible with zer within the errors. Certainly, this 
conclusion should be considered as still preliminary, since the results on the sea contributions 
are obtained for the restricted Bjorken x region. Nevertheless, its degree of reliability is high 
enough due to the discussed above advantage of approximations ( |25|) . (|27|) to the full moments, 
which become especially good for the wide COMPASS x region we deal with. Having in mind 
the surprisingly small [19] [11] values of AiG and Ais it seems that we now became still more 
close to the minimal "retro" picture of the proton spin puzzle, where only helicity PDFs and 
orbital moments of valence quarks compose the nucleon spin. 

Now we are waiting for the new COMPASS data with the planed 180 GeV muon beam, 
which should allow to reach still smaller x values (we expect for the bottom boundary about 
a = 0.003 instead of a = 0.004) and, thereby, increase the reliability of the presented results. 

The authors are grateful to N. Akopov, A. Efremov, O. Ivanov, A. Korzenev, A. Kotikov, 
V. Krivokhizhin, A. Nagaytsev, A. Olshevsky, G. Piragino, G. Pontecorvo, I. Savin, A. Sidorov, 
O. Teryaev, R. Windmolders, and E. Zemlyanichkina for fruitful discussions. 



Appendix 

Combined analysis of data coming from experiments with different binnings 

For brevity and readability we clarify our procedure of combined analysis on a simple ex- 
ample. 

9 While the values of Aid listed in in Tables [SJ El are just zeros within the errors, the values of Aiu are 
compatible with zero within 2a. 

10 Certainly, because of evolution these quantities still can deviate from zero at values of Q 2 distinct from 
considered here Q 2 = 10 GeV 2 . However, in the wide range of Q 2 really available to experiment they are still 
negligible. 

11 Notice that similar conclusion was made in the recent COMPASS paper [IB], where the truncated moment 
of sea u quark is zero within the errors, while the moment of d quark very slightly differs from zero. However, 
first, the analysis in [16] was performed only in LO QCD. Second, the moments studied there were truncated to 
the rather narrow region 0.004 < x < 0.3 (because of lack of the deuteron data in the last two bins). Third, all 
set of FFs has been used, while we apply only well known pion FFs. Fourth, the corrections due to evolution 
were not taken into account in [TB] (approximation (TT5| was applied). 
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Let us suppose that we would like to combine the data on the DIS virtual photon spin 
asymmetry 

A x ~g x /F 1 , 

coming from two independent experiments (expl and exp2) with different binnings. Of impor- 
tance is that our final destination is the integral quantity - truncated first moment - 

rb rb 

^[ a ,b] = / dxgt= dxA x F x (A.l) 

J a J a 

of the structure function g%, where the unpolarized structure function F\ = F2/2x(l + R) is 
considered as an already knowrrn. continuous function of x; a and b are the lowest and highest 
x boundaries available at least to one of the considered experiments (a = 0.004, b = 0.7 for the 
combined analysis of COMPASS and HERMES data we deal with). Of importance is also that, 
just as in Ref. [7], we apply the advanced approximation of the integral by the sum over bins, 
which reproduces much better the real integral value than the "middle point" approximation 
(see discussion around Eq. (16) in Ref. [7]). For the example considered here, it means that 
we approximate by a constant within the bin only the measured asymmetry A x . 

Ax(x) = bms A x {{xi)) 9(x - Xi-i)9(xi - x), 

where Ax((xi)) is the value of asymmetry in i-th bin, 2 = Q, x N bins = b and 9(x) is the usual 
step function, while the known unpolarized input F x (x) remains a continuous function within 
each bin. Thus, instead of the rather crude approximation ( "middle point" approximation) 
r[ 0i 6] ~ J2i=T s Ai((xi)) Fi((xi)) to Eq. flA.ip we use the improved approximation 

r [a , b] = Y!!l7 M{Xl)) r dxF ^ x )- ( A - 2 ) 

Let us now consider the case where some bin [a, (3] of expl is covered by two bins [a, 7] and 
[7,/?] of exp2. Then, to safely combine the respective data with Eqs. (120|) . (12TT) we apply the 
following procedure. First, we artificially divide bin [a, (3} of expl into two bins [a, 7], [7,/?] 
and assign to each of the two new bins a pseudo-measurement of the asymmetry A x . These 

are A ± 5 A and A ± 5 A for the artificial bins [a, 7] and [7, /3], respectively, while the measured 
asymmetry in the initial real bin [a, (3] of expl is just A ± 8 A. Then, we impose the natural 
requirement on our imaginary measurements: 

A = A = A. (A.3) 

Now the pseudo-errors SA, 5 A should be properly adjusted. We do it by imposing the necessary 
condition that the constructed system of pseudo-measurements must be statistically equivalent 
to the system of real measurements within expl. Thus, a first relation connecting the errors 
SA and (6 A, 5 A) is obtained by applying equation (I2TT) : 

(5A) 2 = ^. (A.4) 

1/(5A) 2 + 1/(5A) 2 



12 It is calculated in NLO QCD by using the known expression for DIS Wilson coefficients, some parametriza- 
tion on unpolarized PDFs and the known parameterization on R. 
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At the same time one can see that due to Eq. (IA.30 the second equation ( 120|) of the statistical 
addition just transforms into identity. 

A second relation for the errors is obtained by taking into account that the requirement of 
statistical equivalence should also be satisfied in the calculation of the integral quantity 
we are interested in. In the example considered here it means that the quantity (entering the 
sum in Eq. (|A.2|1 ) r^jli = Af dx Fi and its error ^Tu^li should not change when the bin 
[a, (3} of expl is divided into two pseudo-bins [a, 7] and [7, (3], from which the same quantity is 
derived: 



(A.5) 



r^ln = A dxF t + A dxF 1 ±5T[ a ^]\ n . 

J a J 7 

The requirement of the central values equality 

r[a,/3]|i = r^jin 

is satisfied automatically by virtue of Eq. (jA.30 . while the requirement 



5T 



5T 



a,/8]|n, 



is written as 



13 \ 2 

dx Fx 



,) A I <!.r F ) = \^5AJ dxF x 
Solving the system (1A.4I) . (1A.8j) one obtains 



6A 



f3 \ 2 

dx Fx 



(A.6) 



(A.7) 



(A.8) 



5A = 5A\ 



SldxF, 



5A = 5A, 



Ja dxF i 



(A.9) 



We operate in the same way every time when some bins of expl and exp2 do not coincide 
with each other and, after that, safely apply Eqs. (120]) . ( 12T]) to combine the data of both 
experiments. 

It is easy to see that the generalization of the proposed procedure is straightforward in the 
case of difference SIDIS asymmetries (or in the case of any other asymmetry we deal with). 
For these asymmetries the role of the integral quantity T in the above procedure is played by 
the integral quantities A^" (see (IT2])). giving direct access to the n-th moments of valence 
PDFs we are interested in (see (fTTJ)- 
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Figure 1: Pion difference asymmetries A 7 ^ 



and Aj 



at Q 2 = lOGeV 2 , constructed with Eq. (ITT)) 



from the COMPASS data on Al and A% in the regions 0.004 < x < 0.7 and 0.004 < x < 0.3, respectively. 




and Aj 



Figure 2: Pion difference asymmetries A™ 
from the HERMES data on A^ d in the region 0.023 < x < 0.6. 



at Q 2 = lOGeV 2 , constructed with Eq. (JTTJ) 
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Figure 3: Pion difference asymmetries A* ~ v and A d ~ w at Q 2 = lOGeV 2 , constructed with Eq. ([TTJl 
from the COMPASS data on A* >d in the regions 0.004 < x < 0.7, 0.004 < x < 0.3 (for proton and deuteron 
targets, respectively), and HERMES data on A* d in the region 0.2 < x < 0.6 (last three bins of HERMES). 
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Figure 4: Pion difference asymmetries A* ~* and A d * at Q 2 = 10 GeV, constructed with Eq. dTTJ) from 

the COMPASS and HERMES data on A£ d in the entire x-regions accessible for measurement (0.004 < x < 0.7, 
0.004 < x < 0.3 for COMPASS and 0.023' < x < 0.6 for HERMES). 
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